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In this paper we study the Macroscopic Quantum Oscillation ( MQO) effect in ferromagnetic 
single domain magnets with a magnetic field applied along the hard anistropy axis. The level 
splitting for the ground state, derived with the conventional instanton method, oscillates with the 
external field and is quenched at some field values. A formula for quantum tunneling at excited 
levels is also obtained. The existence of topological phase accounts for this kind of oscillation and 
the corresponding thermodynamical quantities exhibit similar interference effects which resembles to 
some extent the electron quantum phase interference induced by gauge potential in the Aharonov- 
Bohm effect and the O-vacuum in Yang-Mills field theory.. 
PACS number(s): 75.10.Jm, 03.65.Sq, 73.40.Gk, 75.30.Gw. 



O 

Os 
OS 



I. INTRODUCTION 



The remarkable macroscopic quantum phenomena ( MQP ) of nano-magnetic particles at sufficiently low temper- 
atures stimulated considerable research activities in recent years jj] In single-domain particles, MQP represents 
quantum transitions of the magnetic moment M between easy directions in a single-domain ferromagnetic (FM) grain 
and of the Neel vector L in a quantum superposition of AFM sublattices. It has been pointed out that there is a 
crossover at Tg, above which the transition process is dominated by thermal hoping and the transition rate follows 
Arrhenius law r = vq exp[— U /fcsT] where U is the energy barrier and vq is an attempt frequency, below Tb the tran- 
sitions occur by quantum tunneling and a temperature independent rate T = z^oexp[— B] is expected. The crossover 
temperature is defined as Tb — "§7^7- Macroscopic quantum coherence (MQC) phenomena have been studied exten- 
— i \ sively for its exotic characters far from that of classical systems. In 1992, Awschalom et al. M observed evidences of 
resonant tunneling in dense ferritin solutions in which each particle contains about 4500 F(? + spin 5/2 ions below 
200mK. Recently remarkable quantization phenomena were observed in MriyiAc J5| that are also understood as a 
. kind of resonant magnetization tunneling. And resonant tunneling evidences are also shown in Fes &s well as CtRq 
j — . ' where R stand for Mn or Ni ||[7]]. 

t-H , It has been shown that for a wide range of magnetic systems, quantum tunneling is completely suppressed if the total 
spin is a half integer but is allowed for integer-spin particles which is called the "spin-parity effect" fsHlOf and related 

to Kramers' degeneracy.. The "spin-parity effect" is a direct result of the topological phase e % ^ s |ll| and as such does 
not occur in previous Lagrangian formulations of " micromagnetics" , which study the statics and thermodynamics of 
continuous magnetization fields. In this paper we study a kind of MQC phenomena resulting from the topological 
phase, which we call macroscopic quantum oscillation(MQO). This kind of oscillation has been proposed in both 
ferromagnetic (FM) grains and antiferromagnetic molecular magnets in the presence of magnetic fields [^2|-pd|. It 
1 ■ has been pointed that in ferromagnetic particles with biaxial symmetry and an external field along the hard axis the 
tunneling rate-quenching need not be related to Kramers' degeneracy [O). Using spin-coherent-state path integrals, 

O ' 2Kl \ / 1_ ^T 

Garg first showed that the tunneling splitting is quenched |L2| whenever H = (S — n — 1/2) where K\, 

K 2 are the anisotropy constants with K\ > K?, S is the total spins of the particles and n is an integer. In this paper 

. £h \ we solve the paths from the classical action and derive a formula for the level splitting as a function of field, from 

which his specific field values can be read off directly. Following ref. ( Q]), we calculate the tunneling dynamics 

of the magnetic moment of FM particles which may be measured via the static magnetization and the specific heat. 

Here we shall study a spin model with biaxial anisotropy in detail and point out that quantum oscillations are of the 

Aharonov-Bohm-type, similar to that of the 9 vacuum in a Yang-Mills field but different from that in Mn^Ac. 

In Section 2 we show that there is Aharonov-Bohm-type MQO in spin systems with a uni-axial anisotropy. In 

section 3 we discuss a spin model with bi-axial anisotropy and find that the energy level splitting of the ground 

state varies with a vacuum angle 0, which is a good realization of the O vacuum in condensed matter. The instanton 

solution is obtained through the standard procedure and the energy splitting agrees excellently with numerical results. 

We derive a compact formula for the level-splitting induced by tunneling which is valid for the region of low energy 

levels. The results of tunneling effects association with the excited states based on the LSZ (Lehmann, Symanzik 
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and Zimmermann) procedure of field theory are restricted in the low energy region H| [if]. Finally we study the 
thermodynamic properties of the tunneling effect in Section 5. 

II. MQO IN MAGNETIC GRAINS WITH A EASY-AXIS ANISOTROPY 

We consider the following Hamiltonian operator of a ferromagnetic particle with uniaxis anisotropy [ fL9| 

H = Ks 2 z - g f i B Hs z (1) 

which describes uniaxis anisotropy and an magnetic field along hard axis with the anisotropy constants K > 0. The 
spin operators si, i — x,y, z, obey the usual commutation relation [si,Sj-] = ieijk'sk (using natural units throughout). 
The Hamiltonian is exactly diagonal in terms of the eigenstates s z and we have E m — Km 2 — g^sHm. Instead of 
using the result, we express the partition function as a spin-coherent-state path integral for large spins S> 1 |20| 

Z = T r e~ m = J V{^i(n)}e- SE (2) 

where V{fi(n)} = UtL'i 1 sinf>k ^ kMk and 

S E = dT[iS<j>(l-cos9) + KS(S+l)cos 2 6-gfi B HScose}. (3) 
Jo 

Although the method from spin-coherent-state seems redundant, our purpose in this part are to compare this simple 
model with the model in next part. After integrating over cos#, we map the spin system onto a particle problem with 
Lagrangian 

£ = ^+<e* (4) 

where m = and 9 = S ^1 — 9 £g g ^ . This Lagrangian is just similar to that of eq.(3.1) in ref.( |2l|| ). The last term 
of C is the total imaginary time derivative and has no effect on the classical equation of motion, though it alters the 
canonical momentum II^ = m (f> -HO. This magnetic model with uniaxial anisotropy in an external field resembles an 
electron moving on a conducting ring that is crossed by a flux $ which shows remarkable Aharonov-Bohm effect [pl| . 
With imaginary time r = it variable, the periodic instanton configurations are classical solutions which minimize the 
Euclidean action under the boundary condition <f> n (r + (3) — 4>n{T) + 27m 

<Pn = -jf T > ( 5 ) 
Se = son 2 + i2?i7rO 

where n is the winding number characterizing homotopically nonequivalent class and sq = i\ ksT '/ 'K. The Euclidean 
functional integral of the partition function contains an additional summation over the homotopic number 

oo 

Z= Zn = e 3 (7re,e X p(-s )) (6) 

n— — oo 

where <d3(v,q) is the Jacobi theta function oscillating with 8. In the following we shall be interested only in that 
part of the ground state energy which depends on the Q 

E = -k B T\nZ as T ->■ (7) 

\2 



4a 2m 

where {{a;}} is the fractional part of x to the nearest integer. This simple model shows remarkable similarity to 9 
vacuum of the non-Abelian gauge field - both of them are a manifestation of vacuum-vacuum tunneling and angle 
can be made as an external variable parameter in stead of a dynamic one. 
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^From the partition function we find the oscillation of the magnetic momentum with external field H z along 
the hard axis 

dM k B Td 2 \nZ 



OH dH 2 



r l 1 



-.2 



2 + 2^ 



{6} 



as k B Tm <C 1 



-> ^ miB ) - 2k B T (HEtE.) 2 cos (27r9) exp(-s /?i) as k B Tm > 1. 
2K V K J 

where {x} is the integer part of x to the nearest integer. The result shows the regular intervals of magnetic field 
AH = j^. Fig.l shows the magnetization curves at zero temperature M — H : curve I is magnetization curve for 
magnetic particles with integer spin, curve II with half integer spin and curve III is the classical one which has no 
oscillations. This kind of quantum behaviors is essentially a kind of " Aharonov-Bohm effect" in mesoscopic magnetic 
systems resulting from Berry phase. Recently quantum steps in the hysteresis loop of M n \ 2 Ac crystal have been 
observed at regular intervals of magnetic field with period AH = -j^- However the Hamiltonian of M n \ 2 Ac is 
H = —K's 2 — g^j, B H's z + H' where H' is the perturbation term and can't communicate with other terms, such as s± 
and spin-phonon interaction. Even without the term H' , this model can't be mapped onto the above model due to the 
opposite sign of mass, m = ^ . The difference between the two kinds of MQPs is obvious : in the model here MQP 
is essential as a kind of "Aharonov-Bohm effect" in mesoscopic magnetic systems from Berry phase with a period 
AH = jjj^', while in M n \ 2 Ac MQC is a effect of energy level-crossing and the period AH = -^-^ ( The perturbation 
term H' is the key to understand the exotic phenomenon in MnizAc). 

III. MQO IN MAGNETIC GRAINS WITH A BIAXIAL ANISOTROPY 

Let us consider the following Hamiltonian operator of a ferromagnetic particle with XOY easy plane anisotropy and 
x easy axis with a magnetic field applied along z axis fl2]| 

H = K x sl + K 2 s 2 y - gf i B Hs z (9) 
where K\ > K 2 > 0. Starting from the coherent state representation of the time evolution operator with Hamiltonian 

E(e, 4>) = K x cos 2 6 + K 2 sin 2 6 sin 2 - g f i B H cos 6 (10) 
and integrating cos#, we map the spin system onto a particle problem with the effective Lagrangian 

£*tt = \™{<t>) $ -v{<j>) + 6(0) (11) 

where the position dependent mass m(<p) and angle 0(0) are defined by 

771(0) = 5 , (12) 

VW 2tfi(l- Asin 2 0)' V ' 

6(0) =S(1- H 2 ), 

2AiA(l — Asm 0) 

wtx\ v qiq^ t\ ■ 2 A (9VbH) 2 Asin 2 

1/(0) = K 2 S [S + 1) sin - — — 

4iii(l — Asm 0) 

respectively. Here S is the total spin of the particles and A = Jp. The term 6(0) has no effect on the classical 
motion equation and comes from the Berry phase. 

A. Level Splitting of ground states 

Now we derive the level splitting of a giant spin particle in the large spin limit. Instantons in field theory of + 1 
dimensions viewed as pseudoparticles with trajectories existing in barriers, are therefore responsible for tunneling. 
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Since instanton trajectories interpolate between degenerate vacua and satisfy vacuum boundary conditions, the in- 
stanton method is suitable only for the calculation of tunneling parameter Aeo between neighboring vacua. In the 
following we first consider tunneling at the vacuum level (i. e. m = 0), which leads to the level splitting of the ground 
state energy, i.e. 2 Aeo. Passing to imaginary time by Wick rotation r = it, 3 = >T, the amplitude for tunneling from 
the initial well, say that with n = (and 0, = 0), to the neighboring well with n = 1 (and 4>f = vr), is obtained by 



< 0, $„ = 7T 



-13 ft 



0, $„+i = >= K E {4>f = vr, p\ & = 0, 0) 



(13) 



where | 0, > is the perturbation wave function of the ground state at the well denoted by $ n . Since the classical 

■ 2 

solution (f> c satisfies the equation of motion ^m((f> c ) </> — V(0 C ) = 0, we have 



arccos 




(14) 



where uj = y/AK^S {S + 1) (1 - u) and u = 4 ^s(s+i) with 

non-vanishing topological charge and the Euclidean 

action is 



(15) 



where ^ 



(l-A)(l-«) 



■y/S{S+l) 



and u < (1 — A). When H — > 0, we have 



1 — tanh ujqT 
1 — A tanh 2 loqt 



(16) 



which is exactly the vacuum instanton solution derived previously in the absence of magnetic field P, |17||lq ] . 

The prefactor can be evaluated with the stationary phase method by expanding (j) around the instanton trajectory 
4> c such that <j) = 4> c + 77, where 77 is the small fluctuation with boundary conditions tj(/3) = 77(0) = 0. Up to the 
one-loop approximation we have 



where 



T>r\e 



-SS E 



l V (0)=0 

is the fluctuation functional integral with the fluctuation action 



SSe = 



rjMrjdT 







where 



with 



M 



1 d (A\ d 



(17) 



(18) 



(19) 



(20) 



v(4>o) 



J4> c - -m"(0 c )^ + V"{4> c )\ 



dm(<j>) 



(21) 

As in the usual method of 



Here M(0 C ) is the operator of the second variation of the action and ml {(f), 

evaluating the fluctuation integral /, we expand the fluctuation variable r\ in terms of the eigenmodes of M and set 
i] = Y, n C n 4>n, where ip n denotes the n-th eigenfunction of M, and express the result of the integration as an inverse 
square root of the determinant of M. In view of the time translation symmetry of the equation of motion, the functional 
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integral ICe is not well defined when expanded about the classical solutions. The translational symmetry results in 
zero eigenmodes of the second variation operator M of the action. This problem can be cured by the Faddeev-Popov 



procedure J23| . Following the procedure of refs. ]24 . the one instanton contribution to the propagator in the one-loop 
approximation is calculated as 

JC£ 5 =/3-(l-u- (\ - u)i S 2 K 2 e- ulo(s/2 e- s ce-™ B (22) 

To obtain the desired result proportional to sinh(/3 A Eq/2) in cq.([T^), the contributions of the infinite number 
of instanton and anti- instanton pairs to the one instanton contribution have to be taken into account. Interactions 
among instantons and anti-instantons are neglected in the dilute instanton-gas approximation. Summing over all 
contributions of instantons and anti-instantons, the final result of the propagator is found to be 

£b = A*(-)* e~^°/ 2 (23) 

7T 

m+ro odd -. 

x (2D) m+n e -Sc(m+n) e -me(m-n) 

^ mini 

n,m>0 

= A^-)5(1- W )^e"^°/ 2 

7T 



sinh[/3Z? cos (ttQ) e 



where D = 2 3 { (1 ^^ )7r (1 - u)i \iS§ and 9 = S (l - 2S ^" x) i/^ j ■ Compared with esq.© the level shift is seen 
to be 2Aeo cos (-7T0) (note that this is the shift of a single level), the tunneling splitting is 

AE = 2Ae cos (ttG) (24) 

where 

Aeo = 8{ - KlK \. }*(1 u)iA*Sie- s ^. (25) 
(1 — u — \)n 

The tunneling splitting is quenched (A£ = 0) whenever G = k + 1/2 or H = 2y ^ gl (S ~ k - 1/2) /S where k 
is an integer and this is the result of topological phase : there is destructive interference between any path 4>{t) 
and Q(t) with its time reversal pair —4>{t) and 7r — 9(t). The result is in good agreement with earlier works by 
Garg [JL2J. The original Hamiltonian operator in a magnetic field along z axis has no time- reversal symmetry for 
g^BH's z — > — g/is-ffsz as t — > — t where t is time. For this reason, this quenching in magnetic field has nothing to do 
with Kramers' degeneracy and is different from the quenching with half-integer spins in the absence of magnetic field. 



B. Level splitting of low-lying levels and the LSZ method 

A formula suitable for a quantitative of the tunneling effect at excited levels is obtained with the Lehmann, Symanzik 
and Zimmermann (LSZ) method in field theory in which the tunneling is viewed as the transition of n bosons induced 
by the usual(vacuum) instanton. The idea of a tunneling transition from one side of a potential barrier to the other 
has recently also been linked with the LSZ reduction mechanism of a transition from the asymptotic in-state to 
asymptotic out-states. In the following we use the LSZ reduction procedure in a modified way in order to calculate 
the tunneling in the one-instanton sector for the effective potential of the spin system, including the contribution of 
quantum fluctuations up to the one-loop approximation. 

We recall first the case of a one-dimensional harmonic oscillator described by the Hamiltonian 

H = \p 2 + \^lq 2 (26) 

for mass m = 1. Here q and p are dynamical observable which become operators when subjected to the Heisenberg 
algebra of ordinary canonical quantization. The solution of the Heisenberg equation of motion, q + uJ^q = 0, then 
becomes 

q{t) = — =[ae- lUat + a) e luJot ] (27) 
V2^o 
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where a, are time-independent operators defined by the initial (t = 0) values of g and p, i.e. 

q(0) = -jj=[a + a% (28) 

V2wo 

The operators a, a) can be obtained from q(t),p(t) = q(t). Thus 

at = |=e- iwot [g(t) + iwo«(t)] (29) 

V2w 



v 7 ^ & 

and a follows with complex conjugation. One should note the extra — in the definition of the symbol when acting 
to the left. Operators of this type are well-known in the literature. Because the "harmonic oscillator" approximation 
is only useful of low lying levels, the results from LSZ are restricted to low lying levels. 

We now consider the (0+1) dimensional theory defined by the Euclidean Lagrangian eq.(|l2j). Crucial aspects of the 
LSZ procedure are its asymptotic conditions which require the theory to have an interpretation in terms of observables 
for stationary in coming and outgoing states. We can simulate such a situation here artificially by imagining the central 
barrier of the potential to be extremely high and the neighboring wells " — " and " + " on either side to be extremely 
far apart. We therefore construct appropriate functions <f>± (r) which become 

4> + = tt - 4> c , 0_ = 4> c (30) 

such that the interaction fields vanish in their respective asymptotic regions, i.e. 

lim 0+(r)=O, lim 0_(r)=O (31) 

r — >+oo r — > — oo 

The subscripts " — " and " + " here denote the wells with minima at $o = and $i — tt respectively. The Euclidean 
creation and annihilation operators a± and a± which create and annihilate an effective boson in wells " + " and " — " 
respectively are related to the interaction field operators <f>± by 



al(r): = ^e-^ ^- <f> in. (82) 



wo dr 



a ± {T): = -J 2 -^Le^d MT) 
V ujo or 

where := represents definition and tuq — 1/2K\. From the viewpoint of the LSZ method the transition amplitude 
between m-th degenerate eigenstates in any two neighboring wells (here for n = 0, 1) is viewed as the transition of m 
bosons induced by the instanton of eq.dlj) and is related to the tunneling parameter Ae m by 



Af A =< m,$ 1 | e - /Mr |m,*o >= e~ 0em sinh/3Ae m (33) 

where | m, $„ > is the perturbation wave function of the excited level at the well denoted by <& n The transition 
amplitude as well as the ^-matrix can be related to the Green's function through the procedure known as the LSZ 
reduction technique. To this end we rewrite the transition amplitude as 

Af ti = S^e- " (34) 

with S'-matrix element 

S%= lim ^(0\a + ( T i)...a + (T{)aUrl)...aUr l m )\0) (35) 
J ' t»-»-oo m! 

The S'-matrix element can be evaluated in terms of the Green's function G which arises in its evaluation. Thus 



G 



1 m / 

x/^oo 1=1 \ 



wo <9t/ 




9r? 



G 



(36) 



d 2 G 



( dG dG _ A 



where the 2m.-point Green's function is defined as usual, i.e. 

G = (0\$+(t£) . . . $ + (r{)^(rt) . . . 0_«)|O) 
We evaluate G by inserting complete sets of states of final and initial field configurations 4>f,<fii- Thus 

G(Tf,T*)=MT f )^-(T l )A% 



(37) 



(38) 



which vanishes in the limit r l — > —oo,t^ — > oo, due to eq. (|3l|). Thus in eq. ( |36| ) the only nonvanishing contribution 
in these limits results from the second derivative, then 



d 2 G 



6 + (r/)fl0_(r«) AO 



07 



dr-f dr l dr? 
the S'-matrix element for the transition of m bosons is thus 

~# + (r/)#_(r/) 



S 




2mo\ 
w J 



drl 



A 



(39) 



(40) 



The transition amplitude between degenerate ground state can be calculated from the definition, eq. (26), in terms 
of the tunneling parameter Aeo,i.e. A®^ — /3Aeo cos (07r) e~P U0 ' 2 . Then we have 



A f,i - ™|9 e 



m —u)Qm 



where q = 



_ 8sA 1/2 (l-u) 3 
1-A-m 



The tunneling parameter at m-th excited state is seen to be 
Ae m = — :? m Ae cos (Gtt) 



(41) 



(42) 



where Ae m is the usual overlap integral or simply the level shift due to tunneling through any one of the barriers. 

ssA^q-u) 3 / 2 



For q 



l-A- 



(gfj, B H)- 
iK, 



3> 1 the tunneling effect is much more obvious of higher excited states. The energy levels are 



£Ae m cos (On) 



(43) 



+ (m + h)oJo and ^ is an integer and here can take only either of the two values "0" and "T 



where e T , 

In the context of these investigations the usual terminology of MQC refers to the resonance between neighboring 
degenerate wells. There is an essential difficulty related to the existing theory of quantum tunneling itself in the 
absence of an external magnetic field. The argument of the WKB exponential of the tunneling for a ferromagnetic 
particle is 2\/Xs with A = K\ and Ki being the hard and medium axis energies. Recently, evidence for resonant 
tunneling were observed in Mn^Ac, Fes and CtRq ( R equals Mn or Ni) (5|-0|, in which magnetic clusters have 
very small spins (S = 10 — 20) . We show some data for a ferromagnetic particle in order to demonstrate the tunneling 
spliting in magnetic fields H with respect to the result of tunneling for half integer-spin particles in Fig. 2(a) and half 
integer one in Fig. 2(b). The result of the instanton method is plotted with the solid line and the numerical result 
from exact diagonalization is shown by the dotted line. The level splitting of low-lying levels can be also obtained 
with periodic instanton method. For the comparison see the appendices I. 



IV. OBSERVATION OF MQO 

In this part we discuss the possible relevance to experimental test for the topological phase interference or spin-parity 
effect in the single domain FM particles. 
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The direct way to observe MQOs is to observe the resonant tunneling frequency. Let's denote | n >+ and | n >_ be 
the eigenstates of the same energy E n in the (separate) right and left wells, respectively. Due to quantum tunneling, 
the degeneracy is lifted and there emerges an energy level splitting AE n . The new eigenstates are | n > — ^=(| n > + 

— | n >_) with odd parity and | n > e = n > + + \ n >_) with even parity. If we choose | ^(t = 0) >+ = | n > + as 
initial state, the probability of the system being in the other well is | Ck n (t) | 2 =| - < n(t) \ n > + | 2 = sin 2 [ A "^" f ]. The 
resonant tunneling frequencies V — ^j^- — ^ | cos (ttQ) \ Ae q n now with magnetic field. 

Since the topological phase interference can be reflected in thermodynamic quantities, it is reasonable to study the 
thermodynamic properties ( such as the specific heat and the magnetic susceptibility ). Then we discuss thermody- 
namic behavior of this model at very low temperature T ~ Aeo/fcs. At such low temperature the system can be 
regarded as a typical two-level system 



m=0,£ 



-pea 



e ~ cos(G7r) Aeo _|_ g/3 cos(©7r)Aeo 



(44) 



^From the partition function we can calculate the special heat as C„ = ks 



2 cos(e7r)Ae 
fe^T 



sech 2 



2 cos(e7r)A£ 
k7T 



which 



exhibits a characteristic Schottky anomaly |14| shown in Fig. 3. The height of the peak is a constant 0.64fcs and the 
correspond temperature T max exhibits oscillations with a period 2Kl S^ 1 ~ A 



1.2 

T max = — | cos(67r) | Ae . 
Kb 



(45) 



The magnetization M and its derivative with respect to magnetic field are found to be 



M 



(gUB ) H 
2K X 



x tanh 



2 cos (07r) ^ - 2Ae sin (Gtt) J| 
2cos(67r) Ae \ 



(46) 



(gUB) H sin(297r) A 2 e g/iB^ 



2K X 



and 



DM _ {gii B ? _ 2cos(2e7r) A 2 e 
~dR ~ 2K 1 k^f 



2K 1 SVT^X 



(47) 



respectively. The period of oscillation in |^ is 2Kl g ^ X ■ 

When the temperature is higher Aeo ^ ksT < loq, the exited energy levels may give contribution to the partition 
function and we have the partition function as the following form 



Z = ^exp{-/3 



+{m + tt)^o + 2Ae m cos(97r + ^tt) 
4K i Z 



(48) 



Z + 2 [fi Ae cos (6tt)] exp 



(g^sHY 



2qe 



-/3u>/2 



where Iq [x] = 



is the modified Bessel function. We define a characteristic temperature f that is solution 



to the following equation qe u °/ 2J;bT ~ 1. Only below T ■ 
^From the partition function we obtain the free energy 

1 

1 
1 

1 



2k B ln(g) 



can the system regarded as a two-level system. 



F 



\nZ 



(49) 



2[Ae cos(87r)] 2 F(T, H) 
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where F (T, H) = — —, jr-^p — 4- . The magnetization M and its derivative with respect to magnetic field are found 

to be 

(g/XB) 2 g 2^^(29^% 

2ifi 2iY 1 S , % /T^A V ; V ; 

and 

dM = / gfigir \ 2 

d# 2ifi v y V2Xi5 N /T^Ay' v ; V ' 

respectively. 

Our method is only useful at low temperatures and it is not reliable above the crossover temperature Tb ■ Roughly 
speaking the picture of crossover phenomena of transition processes through the barrier can be outlined as follows. 
At very low temperatures (near zero) the transition is dominated by instantons or bounces which satisfy the vacuum 
boundary condition and are responsible for tunneling at ground state. At high temperature the transition is mainly 
due to the thermal activation (over barrier transition) and follows the Arrehenius law(~ exp[— ]) with AJ7 being 
the barrier height. The transition at high temperature is induced by sphelaron from the viewpoint of pseudoparticle 
method. The crossover occurs at certain temperature Tb above which there is only sphelaron and below which only 
periodic instantons contribute to the imaginary part of partition functions. Above the crossover temperature the 
effects of the energy levels above the barrier are more important. 

V. CONCLUSION 

We have discussed a spin model with biaxial anisotropy and find the energy splitting of ground state varies with 
an angle that is a good realization of the 9 vacuum structure in condensed matter. In this paper we show a kind of 
macroscopic quantum oscillation ( MQC) phenomena in magnetic grains that is also the direct result of the topological 
phase and point out that some thermodynamic parameters (such as special heat C v or dM/dH) may show oscillations 
in a magnetic particles with a easy-plane anisotropy, occurring in a external field along the symmetry-axis with the 

period AH = . While in a magnetic field along the medium axis, the small particles also exhibit quantum 

interference effects but there is no oscillation of tunneling splitting with the magnetic field. The energy levels splitting 
are read as the following form 

AE = AE A + AE B cos (2Sn) (52) 

where AE& and AEb have no oscillations with magnetic field 

Appendix 1: Level splitting of low- lying levels and periodic instanton method 

It is by now known that the periodic instanton method has became a powerful tool for evaluation of quantum 
tunneling at entire region of energy. The classical solution <f> c satisfies the equation 

\m{4> c ) £ -V(4> c ) = -E cl (53) 

which is obtained from by integration once and the integration constant E c \ > may be viewed as the classical energy 
of the pseudoparticle configuration 



1 — k 2 sn 2 lot 

1 ^2—2 ( 54 ) 

1 — Afe sn z lot 



where the parameters are 



lj \/1 - n 1 2 nl, (55) 

(56) 








and 



n\ = l/[ 



2A 



A^IP 



(57) 



^2 = A/[— + V T + 



2A 



A 4A J 



o = - (u - 4g J 2 Ag - 1) and b = satisfy the condition n\ > 1 > n\ > 0. 

The amplitude for the transition from the left-hand well to the right-hand well at the energy E c i due to instanton 
tunnelling can be written as 



Af$ = (E cl ,3> n+ i \e~ 2HT \ E cl ,$ n ) = 



-2E cl T 



smh(2TAE c i) 



and can also be evaluated with the help of the path-integral method as 

A U = / i J *E al ,n+l( ( t ) f)^E l:l ,n( C t ) i) IC E((f>f,T f ;(f> i ,T l )d(l)fd(l)i 



(58) 



(59) 



where ipEdinifa)) V , s o i,n+i( ( A/0 are wave functions in the n-th, (n + l)-th well respectively which extend into the 
domain of the barrier. AE c i denotes the energy shift due to tunneling and is the key quantity to be calculated. 
Following the procedure of ref. [p~5| — P~S|| , up to the one- loop approximation we have 



K-E = 



■,-Sc 



and the fluctuation integral can be calculated as 

1 



N(n)N(T f ) 



dr 



N 2 {T)m (</>c0~)) 



(60) 



(61) 



where N(t) = ^ is the zero-mode of the periodic instanton solution and the Euclidean action evaluated for the 
periodic instanton trajectory is given by 



S c (ct>(T),<p(-T),T) 



-T 



^m{cj> c ) £ +V ) (h 



where 



W{4>{T)^{-T),E d ) 



W(<t>{T),<l)(-T),T) + 2E cl T + i6 



[K(k) - (1 - Afc 2 )fI(A£; 2 ,fc)] 



(62) 



(63) 



and the limit <fi(—T) — * ax, 4>{T) — > a 2 is understood. Here K(k),H(\k 2 , k) are the complete elliptic integral of the 
first kind and the third kind, respectively. To complete the end-point integration, we choose the wave functions in 
the barrier as the usual WKB form 



V>£ c! ,n+1 (</>/) 



and the amplitude is written 



A E c i 



C 2 



Cexp C-f*' m{<t>) 4> dcj> 
Cexp (- f£ 2 f m(4>) (/> d<p^j Ce -n^ f ) 

- [** m{<t>)<t>d<t> - r 2 m( 

e Ja i e *'/ 



(64) 



N 2 (r f ) 



d,T 



(65) 



_ T N^r)m(Mr)) 
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The end point integration dfa can be evaluated as 



exp 



m((f>) <f> d<f> ] . 

/ <j> {n) 



2T 



(66) 



under the limit — > oi. We expand and = m(0) d</> as power series of [<pf — (T)] with the limit 0/ — > a 2 
up to the second order for the Gaussian approximation, i.e. 



and 



The transition amplitude then is 



S c = S c {<f>{T),</>{-T),T) 



1 <9 2 



(67) 



(68) 



N' J {T)m((p a (T)) 



(69) 



W -2E cl T * N2 ^f >jZ T J^7 } 



-(4> f -4>(T)) 2 



xe e ~" 0,i e 

The renormalization constant C is defined and evaluated as 

1/2 



(T)m(*) 



c = 



r a 2 d(p 



4K"(fc') 



(70) 



where fc 2 = 1 — ft. Summing over contributions from one instanton plus n instanton-anti-instanton pair, the total 
amplitude is given by 



A Li = E A Li +1) = e ~ 2Ec ' T sinh (2TCe- w ) 



n=0 



where the transition amplitude of one instanton plus n pairs is seen to be 

pT PT\ f-T-2 

dn / dr 2 ■ ■ ■ I 

-T J~T J-T 

The generalized level shift formula due to tunneling for given energy E c i is obtained 



4 2 f +1) = J T d Tl p dn--- dr[C} 2n+1 e^ 2n+1 



)W e -2E al T 



AK(k') 



-w 



(71) 



(72) 



(73) 



and from which we have the energy level-spliting 2AE m cos7r9 where 97r = S[l — —2^==]. 
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Figure Captions 

Fig. 1. Magnetization curves at zero temperature for a magnetic particle : curve I (solid line) is magnetization 
curve for magnetic particles with integer spin, curve II (dotted line) with half integer spin and curve III ( dashed line) 
is the classical one. 

Fig. 2 The energy level splitting AE of ground state for a magnetic particle with K\ = 1, A = 0.2 and (a) S = 16; 
(b) S = 16.5. The result of instanton method is plotted with the solid line and the numerical result is shown by the 
dotted line. 

Fig. 3 The Schottky anomaly behavior of two-level system. 
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Fig.1 
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T max =1.2/k B cos(07t)AE(O) 



